Atom chip based generation of entanglement for quantum metrology 
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Atom chips provide a versatile 'quantum lab- 
oratory on a microchip' for experiments with 
ultracold atomic gases [lj. They have been 
used in experiments on diverse topics such as 
low-dimensional quantum gases |2j, cavity quan- 
tum electrodynamics |3j, atom-surface interac- 
tions [4, 5j, and chip-based atomic clocks [6 J and 
interferometers [3 lj. A severe limitation of 
atom chips, however, is that techniques to con- 
trol atomic interactions and to generate entangle- 
ment have not been experimentally available so 
far. Such techniques enable chip-based studies of 
entangled many-body systems and are a key pre- 
requisite for atom chip applications in quantum 
simulations [9|, quantum information processing 
[10] . and quantum metrology Here we report 

experiments where we generate multi-particle en- 
tanglement on an atom chip by controlling elas- 
tic collisional interactions with a state-dependent 
potential [12 \. We employ this technique to gen- 
erate spin-squeezed states of a two-component 
Bose-Einstein condensate [13j and show that they 
are useful for quantum metrology. The observed 
—3.7 db 0.4 dB reduction in spin noise combined 
with the spin coherence imply four-partite en- 
tanglement between the condensate atoms |14j 
and could be used to improve an interferometric 
measurement by —2.5 ±0.6 dB over the standard 
quantum limit [15j. Our data show good agree- 
ment with a dynamical multi-mode simulation 
|16j and allow us to reconstruct the Wigner func- 
tion [17 \ of the spin-squeezed condensate. The 
techniques demonstrated here could be directly 
applied in chip-based atomic clocks which are cur- 
rently being set up [18j. 

In the currently emerging field of quantum metrology 
[TT] , multi-particle entangled states such as spin-squeezed 
states [THH2T] are investigated as a means to improve 
measurement precision beyond the 'standard quantum 
limit' 15j. This limit arises from the quantum noise 
inherent in measurements on a finite number of uncor- 
rected particles and limits today's best atomic clocks 
[22] . Atom chips combine exquisite coherent control with 
a compact and robust setup [23], suggesting their use 



for quantum metrology with portable atomic clocks and 
interferometers. Several techniques to create entangled 
states on atom chips have been proposed [TU l24lf27] , but 
none has been experimentally realized so far. 

The 'one-axis twisting' scheme of [25] in principle al- 
lows to create a huge amount of entanglement in a two- 
component Bose-Einstein condensate (BEC) [T31 [TBI 125] . 
In this scheme, an initially separable (non-entangled) 
state, where each atom is in a superposition of two in- 
ternal states |0) and |1), dynamically evolves into a spin- 
squeezed state in which the condensate atoms are entan- 
gled. This is due to atomic interactions that provide a 
nonlinear term in the Hamiltonian for the BEC inter- 
nal state. In the experiments reported here, we realize 
this scheme on an atom chip. A notable feature is that 
we control the interactions through the wave function 
overlap of the two states [TBJ [3D] in a state-dependent 
microwave potential [12] . This is a new and versatile 
technique for tuning of interactions in a BEC which also 
works in magnetic traps and for atomic state pairs where 
no convenient Feshbach resonances exist. We use such a 
pair, the hyperfine states |0) = \F = l,mj? = —1) and 
|1) = \F — 2,mp = 1) of 87 Rb, which is also employed 
in chip-based atomic clocks with magnetically trapped 
atoms [18]. 

A two-mode model [16] provides a starting point to 
understand how squeezing is created in our experiment. 
The internal state of a BEC of N two-level atoms can be 
described by a collective spin S = X)ili s «' the sum °f 
the individual spins 1/2 of each atom (see Fig. [T^). Its 
component S z = (Ni—No)/2 is half the atom number dif- 
ference between the states and thus directly measurable. 
A 7r/2-pulse applied to a BEC in 10)®^ prepares it in a 
coherent spin state (|0) + \1))® N /2 N / 2 with mean spin 
(S x ) = N/2 and (S y ) = (S z ) = 0. This is a product state 
in which the atoms are uncorrelated and the quantum 
noise is evenly distributed among the spin components 
orthogonal to the mean spin, AS* 2 = A5 2 = A/4, sat- 
isfying the Heisenberg uncertainty relation AS y AS z — 
\(S x )\/2. This noise gives rise to the standard quantum 
limit if the state is used in a Ramsey interferometer such 
as an atomic clock [22 . 

Quantum correlations between the atoms can reduce 
the variance of one spin quadrature in the yz plane at 
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FIG. 1: Spin squeezing and entanglement through controlled interactions on an atom chip, a, Evolution of the BEC 
internal state on the Bloch sphere (6 = for illustration). Starting with all atoms in |0), a 7r/2-pulse prepares a coherent spin 
state with mean spin (S) along x and isotropic quantum noise in the yz-plane (fuzzy red circle). Subsequent nonlinear evolution 
with xSz deforms the noise circle into an ellipse, creating a spin-squeezed state with reduced noise at an angle # m i n . For state 
tomography, a second pulse rotates the state around —x by a variable angle 9, followed by detection of S z - b, Control of the 
nonlinearity \ on the atom chip. \ depends on the difference of intra- and inter-state atomic interactions. Its dependence on 
the normalized density overlap A of the two BEC components is shown, calculated from stationary mode functions in potentials 
of increasing separation, c, Experimental sequence and motion of the two BEC components corresponding to a. In between 
the pulses for internal-state manipulation (green), a state-dependent microwave potential is turned on (blue; pulse durations 
and microwave ramp times exaggerated for clarity). It dynamically splits and recombines the two BEC components, so that 
X > during the time T. The simulated center-of-mass motion of the two states |0) (black) and |1) (red) is shown as a function 
of time. A slightly asymmetric splitting of the potentials results in an asymmetric oscillation. Insets: corresponding BEC 
mode functions </>o and <f>i along the splitting direction at the beginning, in the middle, and at the end of the sequence, d, 
Measured Ramsey fringes in the normalized population difference JV rc i (here, 8 = 7r/2). Between the pulses, 6 is nonzero (see 
Supplementary Information), resulting in spin precession. The splitting and recombination of the BEC modulates the fringe 
contrast. Here, T is scanned beyond the duration of the sequence in a+c, so that up to ~ 2.5 oscillations are induced. The 
simulated contrast C (red) and density overlap A (blue) are shown for comparison. 



the cost of increasing the variance of the orthogonal one, 
resulting in a spin-squeezed state [28] . To quantify its 
usefulness for metrology, one introduces the squeezing 
parameter £ 2 = iVAS 2 ^/(S,) 2 , where AS 2 min 
is the minimal variance of the spin in the yz plane (see 
Fig. [jj,). The normalization by (S x ) 2 takes into account 
that improving interferometric sensitivity requires not 
only reducing noise but also maintaining high interfer- 
ometer contrast C = 2\(S X )\/N. A state with £ 2 < 1 
allows one to overcome the standard quantum limit in 
a Ramsey interferometer by a factor £ with respect to 
the use of an uncorrelated ensemble of atoms [T^j. Fur- 



thermore, £ 2 is an entanglement witness, with £ 2 < 1 
indicating at least bipartite entanglement between the 
condensate atoms |13) . 

We produce spin-squeezed states by means of time evo- 
lution through the 'one-axis twisting' Hamiltonian |28) : 

H/h = 5S z + ns v + x s 2 z , (i) 

which describes our BEC in good approximation |16j . 
The first term in describes spin precession around z 
at the detuning S. The second term describes spin rota- 
tions around an axis S v = (cos ip) S x — (sin tp) S v due to a 
coupling of |0) and |1) with Rabi frequency fl and phase 
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ip. The third, nonlinear term of strength x arises due to 
elastic collisional interactions in the BEC. It 'twists' the 
state on the Bloch sphere (see Fig. [T^,) , resulting in spin 
squeezing and entanglement. 

An essential feature of our experiment is the control of 
this nonlinearity. Its coefficient 

X = ^ (d No Ho + d Nl Hi - d Nl no - d No Hi) {No ), {Nl) (2) 
depends on derivatives of the chemical potentials 




of the two BEC components evaluated at the mean atom 
numbers (N ) = (Ni) = N/2 after the 7r/2-pulse. Here, 
hj is the single-particle Hamiltonian including kinetic en- 
ergy and the trapping potential, and <fij(r) is the spa- 
tial mode function of state \j). The interaction strength 
cjjk = 4irh 2 ajk/m between atoms in \j) and \k) de- 
pends on the corresponding s-wave scattering length a,jk- 
For our states, the three scattering lengths are close, 
a 00 : otoi : an = 100.4 : 97.7 : 95.0. If the two BEC 
modes overlap spatially, <f>i = 4>o, the crossed terms in 
Q with the minus sign compensate the direct terms 
with the plus sign. Thus, by default, x ~ 0. I n or ~ 
der to increase x> we control the overlap of </>o and </>i 
with a state-dependent trapping potential. By spatially 
separating the two modes, the crossed terms c^^o and 
9jv Mi are set to zero and thus x > 0- In Fig. [T]d, x is 
shown as a function of the normalized density overlap 

A = J dr 3 \cj)o\ 2 \<l)i\ 2 / f dr 3 \<j)o\ 4 f dr 3 ]^ 4 , calculated 
from stationary mode functions in traps of increasing sep- 
aration for our experimental parameters (see below). 

Our experimental setup is described in detail in 
(see also Supplementary Information). In short, we use 
an atom chip to prepare pure BECs of TV = 1250 ± 45 
atoms in state |0) in a harmonic magnetic trap with lon- 
gitudinal (axial) trap frequency /i on g = 109 Hz (/ ax = 
500 Hz). We couple |0) and |1) with O/2-k = 2.1 kHz us- 
ing microwave+rf radiation. The trap minima for |0) and 
|1) can be shifted with respect to each other along the 
longitudinal direction using a chip-based state-dependent 
microwave potential. For detection, we use state-selective 
absorption imaging with a carefully calibrated imaging 
system (see Supplementary Information). It allows us to 
detect both Nq and N\ with good accuracy in a single 
experimental run. 

Our experimental sequence for squeezing (Fig. [l^i+c) 
starts with a resonant 7r/2-pulse of duration 120 /is 
to prepare the coherent spin state. During the pulse, 
O 3> X-/V so that the nonlinearity can be neglected. Af- 
ter the pulse, we squeeze the state by turning on x f° r 
a well-defined time by spatially splitting and recombin- 
ing the two components of the BEC in the following 
way (see Fig. fit). The microwave potential is turned 



on within 50 ^s resulting in an abrupt separation of 
the trap minima for |0) and |1) by s = 0.52 fim. The 
two components of the BEC start to perform one oscil- 
lation in their respective potentials. During the oscilla- 
tion, which is strongly influenced by mean-field effects, 
the mode functions <f>$ and 4>\ almost completely sepa- 
rate so that x = 1.5 s _1 at maximum separation. After a 
time T = 12.7 ms the states overlap again, the microwave 
potential is switched off within 50 fis, and the squeezing 
dynamics, as well as the relative atomic motion, stops. 
This value of T nearly coincides with the 'best squeezing 
time' expected from the two-mode model 16 . We ana- 
lyze the produced state by performing state tomography. 
With the mean spin along x, we measure the transverse 
spin components Sg = (cos9)S z — (sin 9) S y along any 
angle 6 by rotating the state vector in the yz-plane by 
that angle prior to detection of S z — (Ni~N )/2. This is 
done by applying a second pulse for a duration tq =8/0, 
and with a phase tp = -k (ip = 0) for turning clockwise 
(counterclockwise) . 

Figure [2^, shows the noise in Sg obtained from a large 
number of such measurements as a function of 9. Data 
for a squeezed state is shown in comparison with data 
for a coherent spin state where the traps were not sepa- 
rated during the sequence (reference measurement). We 
plot the normalized variance A n Sg = AASg/(N), so 
that A n Sg = dB corresponds to the standard quantum 
limit. In the squeezed state, the spin noise A n Sg falls sig- 
nificantly below the standard quantum limit, reaching a 
minimum of AS" 2 = -3.7±0.4 dB at 9 min = 6°. The cor- 
responding interference contrast is C = (88 ± 3) %. This 
results in a squeezing parameter of £ 2 = —2.5 ± 0.6 dB, 
proving that the state is useful for quantum metrology 
and that the condensate atoms are entangled. The refer- 
ence measurement, by contrast, stays above the standard 
quantum limit for all values of 9. 

Entanglement in the BEC has been defined as the non- 
separability of the iV-particle density matrix [13l [M] . An 
intriguing question concerns the depth of entanglement: 
How large must the clusters of entangled atoms be at 
least in order to produce the observed squeezing? In [T3] , 
a procedure to determine the depth of entanglement from 
the measured spin noise reduction and mean spin length 
is described. Our data falls below the spin-3/2-curve in 
Fig. 1 of [2], which is an experimental proof that the 
condensate atoms are entangled in clusters of at least 
4 ± 1 particles. 

During the splitting, internal and motional states of 
the atoms are entangled. This can lead to a decrease 
of C and thus an increase of £ 2 if the recombination is 
not perfect. To find the time of maximum contrast, we 
perform a Ramsey measurement where the second pulse 
area is 9 = n/2 and T is scanned. Figure [TJi shows the 
resulting Ramsey fringes in the normalized population 
difference 7V rcl = (Nt - N a )/N as a function of T. We 
observe a high contrast of C = (88±3) % at T = 12.7 ms, 
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FIG. 2: Spin noise tomography and reconstructed 
Wigner function of the spin-squeezed BEC. a, Ob- 
served spin noise for the spin-squeezed state (solid circles) and 
for a coherent spin state (reference measurement, open cir- 
cles). The normalized variance A n Sg = 4ASg/(N) is shown 
as a function of the turning angle 9 in the yz-plane, with sta- 
tistical error bars. For this graph, we remove photon shot 
noise due to the imaging process as described in the Supple- 
mentary Information. In the squeezed state, a spin-noise re- 
duction of — 3.7±0.4 dB is observed for 9 m m = 6°, correspond- 
ing to £ 2 = — 2.5±0.6 dB of metrologically useful squeezing for 
our Ramsey contrast of C = (88 ±3) %. Solid lines are results 
from our dynamical simulation. Blue: squeezed state with 
losses but without technical noise; red: squeezed state with 
losses and technical noise; black: reference measurement with 
losses and technical noise, b, Wigner function of the spin- 
squeezed BEC reconstructed from our measurements. The 
black contour line indicates where the Wigner function has 
fallen to l/y/e of its maximum. Squeezed and 'anti-squeezed' 
quadratures are clearly visible. For comparison, the circular 
1 / \fe contour of an ideal coherent spin state is shown. The 
area of the contour line is larger than the area of the circle in- 
dicating that the squeezed state is not a minimum uncertainty 
state anymore. 



indicating large spatial overlap and nearly vanishing rela- 
tive motion of the two states. In the squeezing sequence, 
we turn off the microwave potential at this time, preserv- 
ing the large overlap for subsequent measurements. The 
contrast could be further increased using optimal control 
techniques 25J . In comparison with the data, we show C 
and A as obtained from a simulation for our experimental 
parameters. For an accurate description of our system, 
accounting for both the spatial and the spin dynamics, 
we use the dynamical multi-mode theory developed in 
[To] . It neglects initial thermal excitations and reduces to 
the simple two-mode model described above for the case 



of stationary condensates. The only adjustable parame- 
ter in the simulation is the splitting distance s, which is 
not resolved by our imaging system. The resulting value 
s = 0.52 /im is consistent with a simulation of the poten- 
tial. We observe very good agreement with the measure- 
ment, indicating that the simulation correctly accounts 
for the motion of the BEC in the trap. In Fig [Tj;, the 
simulated dynamics of the mode functions <f>i is shown. 

The spin noise reduction obtained from this simulation 
is shown in Fig. [2^, along with the data. We add the effect 
of particle loss (1,2,3-body) as determined from the two- 
mode Hamiltonian ([!]) in a Master Equation approach 
[TBI [29] . as well as several technical noise sources. The 
blue line shows the expected squeezing taking into ac- 
count atom loss but no technical noise. The maximal re- 
duction in variance is —12.8 dB, significantly larger than 
observed. The red line, which describes our data well, 
additionally includes the fluctuations of N, fluctuations 
of the pulse power of 0.5% r.m.s., a fluctuating detun- 
ing of 27T x 40 Hz r.m.s. during the pulses, and phase 
noise of Atp = 8° r.m.s. (see Supplementary Informa- 
tion). All fluctuations are consistent with independent 
measurements. The fluctuating detuning is due to fluc- 
tuating microwave level shifts during the pulses. It is the 
cause for A n Sg > dB at 9 = 180°. The phase noise 
is the main reason why the maximum achieved squeez- 
ing is smaller than the theoretically predicted value. It 
is consistent with technical fluctuations of the magnetic 
trap position in the inhomogeneous microwave near-field. 
Consequently, the phase noise in the reference sequence 
is smaller (Atp = 3°, black line). 

The measured histograms of Se for different angles 9 
are tomographic data that allow us to reconstruct the 
Wigner function W(S y , S z ) of the squeezed BEC [17] us- 
ing the inverse Radon transform (see Supplementary In- 
formation). Figure [2Jd shows the reconstruction. The 
two contour lines indicate where the Wigner functions of 
our squeezed state and of an ideal coherent spin state 
(with the same N and with added imaging noise) have 
fallen to l/\/e of their maximum. The squeezing along 
the direction 6* m j n as well as the 'anti-squeezing' in the 
perpendicular direction can be clearly seen. 

In conclusion, using a novel method to control interac- 
tions with a state-dependent potential, we have demon- 
strated for the first time spin squeezing and multi-particle 
entanglement on an atom chip. We envisage the imple- 
mentation of this technique in portable atomic clocks and 
interferometers operating beyond the standard quantum 
limit. Furthermore, it is a valuable tool for experiments 
on many-body quantum physics and could enable quan- 
tum information processing on atom chips [25] . 

The group around M. Oberthaler has independently 
and simultaneously realized spin squeezing through Fesh- 
bach control of interactions in an optical trap. 



5 



* E-mail: alice.sinatra@lkb.ens.fr 



philipp.treutlein@unibas.ch 



* E-mail: treutlein@lmu.de 
Fortagh, J. & Zimmermann, C. Magnetic microtraps for 
ultracold atoms. Rev. Mod. Phys. 79, 235 (2007). 
Hofferberth, S., Lesanovsky, I., Fischer, B., Schumm, 
T. & Schmiedmayer, J. Non-equilibrium coherence dy- 
namics in one-dimensional Bose gases. Nature 449, 324 

(2007) . 

Colombe, Y. et al. Strong atom-field coupling for Bose- 
Einstein condensates in an optical cavity on a chip. Na- 
ture 450, 272 (2007). 

Lin, Y., Teper, I., Chin, C. & Vuletic, V. Impact of 
the Casimir-Polder potential and Johnson noise on Bose- 
Einstein condensate stability near surfaces. Phys. Rev. 
Lett. 92, 050404 (2004). 

Aigner, S. et al. Long-range order in electronic transport 
through disordered metal films. Science 319, 1226-1229 

(2008) . 

Treutlein, P., Hommelhoff, P., Steinmetz, T., Hansen, 
T. W. & Reichel, J. Coherence in microchip traps. Phys. 
Rev. Lett. 92, 203005 (2004). 

Wang, Y.-J. et al. Atom Michelson interferometer on a 
chip using a Bose-Einstein condensate. Phys. Rev. Lett. 
94, 090405 (2005). 

Schumm, T. et al. Matter-wave interferometry in a dou- 
ble well on an atom chip. Nat. Phys. 1, 57 (2005). 
Lloyd, S. Universal quantum simulators. Science 273, 
1073 (1996). 

DiVincenzo, D. The physical implementation of quantum 
computation. Fortschr. Phys. 48, 771 (2000). 
Giovannetti, V., Lloyd, S. & Maccone, L. Quantum- 
enhanced measurements: Beating the standard quantum 
limit. Science 306, 1330-1336 (2004). 
Bohi, P. et al. Coherent manipulation of Bose-Einstein 
condensates with state-dependent microwave potentials 
on an atom chip. Nat. Phys. 5, 592 (2009). 
S0rensen, A., Duan, L.-M., Cirac, J. I. & Zoller, P. Many- 
particle entanglement with Bose-Einstein condensates. 
Nature 409, 63 (2001). 

S0rensen, A. & M0lmer, K. Entanglement and extreme 
spin squeezing. Phys. Rev. Lett. 86, 4431 (2001). 
Wineland, D. J., Bollinger, J. J., Itano, W. M. & Heinzen, 
D. J. Squeezed atomic states and projection noise in 
spectroscopy. Phys. Rev. A 50, 67-88 (1994). 
Li, Y., Treutlein, P., Reichel, J. & Sinatra, A. Spin 
squeezing in a bimodal condensate: spatial dynamics and 
particle losses. Eur. Phys. J. B 68, 365-381 (2009). 
Wigner, E. On the quantum correction for thermody- 
namic equilibrium. Phys. Rev. 40, 749-759 (1932). 
Rosenbusch, P. Magnetically trapped atoms for compact 
atomic clocks. Appl. Phys. B 95, 227-235 (2009). 
Esteve, J., Gross, C, Weller, A., Giovanazzi, S. & 
Oberthaler, M. K. Squeezing and entanglement in a bose- 
einstein condensate. Nature 455, 1216-1219 (2008). 
Appel, J. et al. Mesoscopic atomic entanglement for pre- 
cision measurements beyond the standard quantum limit. 
Proceedings of the National Academy of Sciences 106, 
10960-10965 (2009). 

Schleier-Smith, M. H., Leroux, I. D. & Vuletic, V. 
Reduced-quantum-uncertainty states of an ensemble of 
two-level atoms. preprint \arXiv. 08 10.2582] (2008). 



[22] Santarelli, G. et al. Quantum projection noise in an 
atomic fountain: A high stability cesium frequency stan- 
dard. Phys. Rev. Lett. 82, 4619-4622 (1999). 

[23] Vogel, A. et al. Bose-Einstein condensates in micrograv- 
ity. Appl. Phys. B 84, 663 (2006). 

[24] Calarco, T. et al. Quantum gates with neutral atoms: 
Controlling collisional interactions in time-dependent 
traps. Phys. Rev. A 61, 022304 (2000). 

[25] Treutlein, P. et al. Microwave potentials and optimal 
control for robust quantum gates on an atom chip. Phys. 
Rev. A 74, 022312 (2006). 

[26] Charron, E., Cirone, M. A., Negretti, A., Schmiedmayer, 
J. & Calarco, T. Theoretical analysis of a realistic atom- 
chip quantum gate. Phys. Rev. A 74, 12308 (2006). 

[27] Zhao, B., Chen, Z., Pan, J. & Schmiedmayer, J. High- 
fidelity entanglement via molecular dissociation in inte- 
grated atom optics. Phys. Rev. A 75, 042312 (2007). 

[28] Kitagawa, M. & Ueda, M. Squeezed spin states. Phys. 
Rev. A 47, 5138-5143 (1993). 

[29] Li, Y., Castin, Y. & Sinatra, A. Optimum spin squeezing 
in Bose-Einstein condensates with particle losses. Phys. 
Rev. Lett. 100, 210401 (2008). 

[30] Poulsen, U. & M0lmer, K. Quantum beam splitter for 
atoms. Phys. Rev. A 65, 33613 (2002). 



Acknowledgements 

It is a pleasure to thank K. M0lmer, J. Reichel, A. 
Smerzi, and A. S0rensen for helpful discussions and J. 
Halimeh for careful reading of the manuscript. This 
work was supported by the Nanosystems Initiative Mu- 
nich. T.W.H. gratefully acknowledges support by the 
Max-Planck- Foundation. 



Author information 

Correspondence and requests for materials should be 
addressed to P.T. (experiment) and A.S. (simulation). 

Author contributions 

A.S. and P.T. jointly conceived the study. M.F.R., 
P.B., and P.T. carried out the experiment and analyzed 
the data. Y.L. and A.S. carried out the simulations. 
All authors discussed the results and contributed to the 
manuscript. 



6 



ATOM CHIP BASED GENERATION OF ENTANGLEMENT FOR QUANTUM METROLOGY: 

SUPPLEMENTARY INFORMATION 



Experimental setup 

The core of our experiment is an atom chip with inte- 
grated microwave guiding structures. It allows us to gen- 
erate state-dependent microwave near-field potentials in 
addition to static magnetic traps. The chip, the prepa- 
ration of Bose- Einstein condensates (BECs), and the use 
of microwave near-field potentials for state-dependent co- 
herent manipulation of ultracold atoms is described in 
detail in pQ. We briefly summarize it in the following, 
highlighting the differences to pQ. 

Our experimental sequence starts by preparing a BEC 
in state |0) without a discernible thermal component in a 
static magnetic trap on the atom chip. Magnetic shield- 
ing and the use of stable current sources ensure stable 
magnetic potentials which allows us to prepare BECs 
with well-defined total atom number N — 1250 ± 45 
through radio-frequency evaporative cooling. The exper- 
iment is performed in a cigar-shaped magnetic trap at a 
distance of 44 /im from the atom chip surface with lon- 
gitudinal (axial) trapping frequency of /i on g = 109 Hz 
(/ax = 500 Hz) and a magnetic field in the trap center of 
B = 3.36 G. 

The BEC internal state is manipulated by coherently 
coupling the two-photon transition |0) O |1) with radio 
frequency and microwave radiation from an off-chip an- 
tenna and horn, respectively. The microwave is tuned 
2ir x 360 kHz above the transition to the intermediate 
state \F — 2,mp — 0), resulting in a two-photon Rabi 
frequency of fi = 2ir x 2.1 kHz. Fig. [3] shows the resulting 
Rabi oscillations for a detuning 5 = from two-photon 
resonance. The efficiency of a 7r-pulse is (96 ± 1) %. 

Figure |4^i shows Ramsey interference fringes between 
|0) and |1) as a function of the delay T between two ir/2- 
pulses. This is a reference measurement taken in a static 
magnetic trap without splitting the condensate, i.e. with 
the BEC in a coherent spin state. The Ramsey contrast 
at T = 12.7 ms is C = (96 ± 1) %. While the pulses are 
applied, the two-photon resonance frequency is shifted by 
A mw = 2tt x 7.6 kHz with respect to the undriven system. 
This is due to differential AC Zeeman level shifts of |0) 
and 1 1) caused by the detuned microwave radiation of the 
two-photon drive [2J. We always adjust the frequency of 
the two-photon drive such that the detuning from two- 
photon resonance is 6 = while the pulse is applied. In 
between the pulses, the phase of the atomic superposition 
state thus evolves at a rate — A mw with respect to the 
two photon drive. This determines the frequency of the 
Ramsey oscillations in Fig. [4ji. 

For state-selective spatial splitting and recombination 
of the two BEC components, we use a microwave near- 
field potential created with an on-chip waveguide. Com- 
pared with pQ, the detuning of the microwave near-field 




preparation time (|js) 

FIG. 3: Rabi oscillations. Resonant Rabi oscillations of the 
relative atom number iV r ei = (Ni — Na)/(Ni + No) recorded 
by varying the duration of the state preparation pulse. The 
efficiency of a 7r-pulse is (96 ± 1) %. The decay with a time 
constant of 15 ms is due to gradients in Q near the structured 
metallic chip surface which imposes boundary conditions on 
the electromagnetic field. 



from the transition \F = 1, mj? =0) O \F = 2, mp = 0) 
is much larger, 2ir x 12 MHz, so that admixtures of other 
states to |0) and |1) are smaller and the states are more 
robust against magnetic field noise. In this configuration, 
both states experience a microwave potential of opposite 
sign, and with different magnitude due to the different 
hyperfine transition strengths. At a microwave power of 
P mw = 120 mW launched into the chip, a splitting of 
the potential minima for the two states of s = 0.52 /xm 
results. 

Figure |4j) shows Ramsey fringes measured in the 
squeezing sequence, i.e. the BEC is split and recombined 
during the time T between the 7r/2-pulses as in Fig. 1 of 
the main text. Turning on the microwave near-field po- 
tential has two effects: The oscillation frequency of the 
Ramsey fringes slightly decreases because of the differ- 
ential energy shift experienced by the two states in the 
potential. More importantly, the fringe contrast is mod- 
ulated by the overlap of the BEC mode functions <j>o and 
(j>\. The contrast at T — 12.7 ms, the time at which the 
squeezed state is analyzed in the experiment, determines 
the length of the mean collective spin: C = 2\(S X )\/N. 
We observe a contrast of C = (88 ± 3) %, smaller than in 
the reference because the overlap of the BEC mode func- 
tions after splitting and recombination is less than unity. 
The difference to the contrast predicted by the dynamical 
simulation of 94 % can most likely be explained by small 
motion in the transverse direction which is excited in the 
experiment but not modeled. The contrast could be in- 
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time between jr/2-pulses T (ms) 

FIG. 4: Ramsey interference fringes. Ramsey interfer- 
ence fringes in the relative atom number i\T re i recorded by 
varying the delay T between two 7r/2-pulses. a, Ramsey 
fringes in the reference sequence in a static magnetic trap. 
The contrast obtained from a sinusoidal fit is C = (96 ± 1) %. 
b, Ramsey fringes in the squeezing sequence with state- 
selective splitting and recombination of the BEC embedded 
between the 7r/2-pulses, as in Fig. Id of the main text. The 
contrast is C = (88 ± 3) %. 

creased by optimal control of the atomic motion rather 
than abrupt switching of the potentials, as described in 



Imaging system 

Compared with [I], we use an improved detection sys- 
tem which allows us to take absorption images of both 
states and determine Nq and N± in a single shot of the 
experiment. We use a back-illuminated deep-depletion 
CCD camera with a quantum efficiency of 90 % at 780 nm 
and fast line transfer. After the squeezing sequence 
described in the main text, the atoms are transferred 
within 30 ms into a relaxed trap with /i on g = 40 Hz and 



/ax = 130 Hz at a distance of 200 from the chip sur- 
face. The trap is switched off and states |1) and |0) are 
imaged after times-of-ffight of 4.6 ms and 6.1 ms, respec- 
tively. State |1) is directly imaged with a <7~-polarized 
resonant laser beam on the F = 2 — > F' = 3 cycling 
transition. After the image is taken, atoms in 1 1) fly out 
of the depth of focus of the imaging system due to the 
photon recoil momentum transferred during the imag- 
ing pulse. Subsequently, state |0) is optically pumped 
into the F = 2 manifold of the ground state using a 
F = 1 —> F' = 2 pumping laser [3] and imaged on the 
F = 2 — s- F' = 3 transition. For both states, the imag- 
ing pulse duration is 40 (is and the imaging intensity is 
/ = 0.8 / sa t, where I sa ± is the saturation intensity on the 
cycling transition. The FWHM diameters of the imaged 
atom clouds are 15 /xm in the vertical and 10 /im in the 
horizontal direction, both larger than the optical reso- 
lution of our imaging system of 4 /im. The maximum 
optical densities in the cloud centers for 600 atoms in 
each state are 1.2 and 1.4, respectively. The good agree- 
ment of the observed Rabi oscillations with the expected 
sinusoidal behavior proves the linearity of our imaging 
system. 

For the correct determination of the fluctuations of Sg 
it is crucial to know the total atom number N = N + Ni 
accurately. We calibrate our imaging system by two in- 
dependent methods. First, following the method of [I], 
we adjust the effective scattering cross section a e s which 
is used to calculate the atom number from the optical 
density of the cloud such that the measured atom num- 
ber is independent of the imaging light intensity. We find 
Ccff = 0.9 <7oi where do is the theoretically expected scat- 
tering cross section on the \F — 2,mp = —2) «-» \F' = 
3, m' F = —3) cycling transition. This is plausible tak- 
ing into account optical pumping in the beginning of the 
imaging pulse and imperfect imaging light polarization. 

An independent test of this calibration can be obtained 
by observing the scaling of projection noise with total 
atom number for a coherent spin state. Figure [5] shows 
the variance AS* 2 measured directly after a 7r/2-pulse as 
a function of the total atom number N. The constant 
offset due to imaging noise is subtracted as described 
below. The observed linear behavior confirms that pro- 
jection noise ASf cx N dominates over technical noise 
which generically scales as AS 12 oc A^ 2 . If we use the 
first method to calibrate N, a fit to the data in Fig. [5] 
with a straight line through the origin yields a slope of 
0.22 ± 0.01. This agrees with the theoretically expected 
slope of 1/4 to better than 15%. The difference lies 
within the error bar of our atom number calibration ac- 
cording to the first method. As the dependence of AS* 2 
on A" can be very accurately determined from Fig. [5j we 
use it to calibrate the total atom number by rescaling A~ 
so that the slope of the linear fit is 1/4. 
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total atom number N 

FIG. 5: Projection noise as a function of atom number. 

A measurement of the variance AS 2 directly after a w/2- pulse 
is shown as a function of total atom number N, where N is 
independently calibrated with the method of i. . The dashed 
line shows the expected linear scaling with a slope of 1/4, 
the blue line is a linear fit to the data which yields a slope 
of 0.22 ± 0.01. The red line is a quadratic fit with AS:? = 
aN + bN 2 . It yields a = 0.22±0.02 and b = (0.06 ±1) x 10 -5 , 
confirming the linear scaling expected for projection noise. 

Data evaluation 

Our experiment runs very stable and we can usually let 
it record data overnight without supervision. The data 
shown in Fig. 2 of the main text is gathered during 7 
measurement nights with between 80 (for 20 < 9 < 90) 
and 370 (for < 9 < 20) experimental realizations per 
point. The error bars are statistical errors assuming a 
normal distribution. 

We analyze the recorded images by counting atoms in 
the two states within two small rectangular regions with a 
typical size of 30 /j,m x 30 /im. This yields Nq and N\ and 
thus S z = (Ni — N )/2. Only shots where the total atom 
number N differs by no more than 150 from the mean 
N = 1250 are used for the analysis and we check that 
a tighter (75 atoms) or wider (250 atoms) post-selection 
does not significantly change the data quality. In the data 
for 90° < 9 < 360° a slow drift oiNi-N is observed. We 
correct for this technical drift by subtracting a filtered 
data set from the respective raw data, using a second 
order Savitzky-Golay filter [5] over 300 shots. 

In Fig. 2a of the main text, we additionally correct AS 2 
for noise in our imaging system. Photon shot noise from 
the imaging beam contributes to the measured fluctua- 
tions of iVo and Ni with a standard deviation of AN 0: p Sn 
and AiVi.psn, respectively. We determine this noise either 
by measuring the apparent atom number fluctuations in 
regions of the image where no atoms are present, by tak- 
ing 'reference shots' where no atoms are prepared in the 
first place, or by calculating the expected shot noise from 



the observed imaging light intensity. All methods yield 
similar results. We correct the variance of S z as 

ASl corr = AS 2 - (ANl psn + A7V 1 2 psn )/4. 

The applied correction due to imaging noise corresponds 
to typically (AJV 2 psii+ A# 1 2 psn ) 1 / 2 /2 « 7 atoms. Finally, 
we normalize the variance to the expected variance for a 
coherent spin state with (No) = (N%) — (N)/2, obtaining 
A n S 2 z = 4 AS 2 co „/(N). Without subtraction of imaging 
noise, we still observe a reduction in the spin fluctuations 
of -2.3 dB. 



Phase noise 

In our experiment, phase noise is the main reason why 
we do not reach the theoretically predicted spin noise 
reduction. It is therefore important to identify and elim- 
inate the main technical sources for this noise. In the 
squeezing sequence with 9 = 6°, where we observe the 
minimum of AS 2 , we measure the dependence of Sg on 
various experimental parameters. From this we can cal- 
culate the sensitivity of the relative phase between the 
atomic state and our two-photon drive on these param- 
eters. We independently determine the technical fluctu- 
ations of these parameters and estimate their contribu- 
tion to the total phase noise of Aip — 8°. Timing jitter 
of our experiment control (sa 100 ps), fluctuations of the 
external magnetic field (reduced by a fi-metal shield to 

0.3 mG), power fluctuations (~ 5 x 10 -3 ) and phase in- 
stabilities (« 0.3°) of the radio frequency and microwave 
generators and amplifiers for the two-photon drive, and 
fluctuations of the on-chip currents together contribute a 
phase noise of about 1°. The microwave power coupled 
into the on-chip waveguide for creating the near-field po- 
tential fluctuates by 60 /jW (5 x 10 -4 ) leading to a phase 
noise of 2° . The remaining phase noise is consistent with 
fluctuations of a current source used to create a magnetic 
field for the static magnetic trap. Its effect on the phase 
is not due to the small differential magnetic moment be- 
tween |0) and |1) but due to a shift of the magnetic trap 
position in the highly inhomogeneous microwave near- 
field which leads to fluctuating microwave level shifts. 

Wigner Function reconstruction 

We reconstruct the Wigner function of the spin- 
squeezed condensate in the following way: For each mea- 
sured 9 £ [—90°, 90°] we create a histogram of Sg and 
fit it with a cubic spline to obtain a smooth curve. We 
then use a filtered back-projection algorithm [6] to per- 
form an inverse Radon transform [7] . The inverse Radon 
transform is derived for classical image reconstruction in 
a plane. In [8], it has been used to reconstruct the Wigner 
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function of squeezed states of the electromagnetic field. 
However, it is generally not suited to reconstruct an arbi- 
trary spin state on the curved Bloch sphere. In our case 
the spin-squeezed state does not 'wrap around' the Bloch 
sphere so that we can locally approximate the Bloch 
sphere by a plane and use the inverse Radon transform. 
We furthermore make a continuum approximation to the 
measured values of S z , which is reasonable as our imag- 
ing system does not have single atom resolution. With 
the experimental method presented here but with a more 
sophisticated analysis [5] , the density matrix of arbitrary 
spin states that spread over the whole Bloch sphere can 
be reconstructed. Quantum state tomography is of inter- 
est as it gives access to measures of entanglement such 
as the quantum Fisher information [TUJ E] , which char- 
acterizes a more general class of states (including states 
with £ 2 > 1) that can be used to overcome the standard 
quantum limit [11 . 
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